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Abstract 

We revisit the model of a system made up of a Brownian quantum oscillator under the influence 
of an external classical force and linearly coupled to an environment made up of many quantum 
oscillators at zero or finite temperature. We show that the HPZ master equation for the reduced 
density matrix derived earlier [B.L. Hu, J. P. Paz, Y. Zhang, Phys. Rev. D 45, 2843 (1992)] with 
coefficients obtained from solutions of integro-differential equations can assume closed functional 
forms for a fairly general class of spectral densities of the environment at arbitrary temperature 
and coupling strength. As an illustration of these new results we solve the corresponding master 
equation and calculate, among other physical quantities, the uncertainty function whose late time 
behavior can be obtained fully. This produces a formula for investigating the standard quantum 
limit which is central to addressing many theoretical issues in macroscopic quantum phenomena 
and experimental concerns related to low temperature precision measurements. We find that any 
initial state always settles down to a Gaussian density matrix whose covariance is determined by 
the thermal reservoir and whose mean is determined by the external force. For more general spectra 
we show that the solution of the master equation can be reduced to solving for the motion of a 
classical parametric oscillator with parametric frequency determined by the unsolved for master 
equation coefficients. States in these systems experience evolution that is parametrically similar 
to the simpler evolution explicitly determined for in the case of Laurent-series spectra. 
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I. INTRODUCTION 

A. Historical Context and New Results 

In this paper we continue the lineage of work on quantum Brownian motion via the 
influence functional path-integral method of Feynman and Vernon[i] used by Caldeira and 
Leggett[ ] to derive a master equation for a high-temperature ohmic environment, which 
corresponds to the Markovian regime. Following this, Caldeira, Cerdeira and Ramaswamy 
(CCR)[ I] derived the Markovian master equation for the system with weak coupling to 
an ohmic bath, that was claimed to be valid at arbitrary temperature (Sec. II D). At the 
same time Unruh and Zurek[ ] derived a more complete and general master equation that 
incorporated a colored noise at finite temperature. Finally Hu, Paz and Zhang (HPZ) [12] 
derived an exact master equation for a general environment (arbitrary temperature and 
spectral density) which preserves the positive definiteness of the reduced density matrix, an 
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important property missed out in many earlier derivations. 

For many calculations, physicists often invoke Markovian master equations as they are 
easier to solve owing to their having compact, functional representations. But one runs 
into trouble if one is interested in low temperatures, short times or non-ohmic baths, which 
likely fall under the non-Markovian regime (see Ref. [21] for a discussion). The HPZ equa- 
tion is capable of dealing with the full range of parameters for a general environment but 
its form is somewhat involved. For example the coefficients are represented by solutions to 
integro-differential equations and multiple integrals over Green's functions. In this paper we 
present the correct (c/. for instance OCR) late-time non-Markovian master equation in com- 
pact functional representation for arbitrary temperature and coupling strength (Sec. II C 2). 
Additionally we obtain the coefficients for a large class of spectral densities and their time 
dependent form at times before the oscillator and reservoir have equilibrated. 

Solutions to the HPZ equation have been attempted before [22] and the role of the master 
equation coefficients (renormalized frequency, dissipation, diffusion, and anomalous diffusion 
coefficients) has been discussed in the past. Here we present complete solutions of the master 
equation. We consider a quantum oscillator under the influence of an external classical force 
that is linearly coupled to a thermal reservoir of quantum oscillators. The particle begins 
decoupled from the reservoir and uncorrelated with it, with a short switch-on time r for 
the coupling. The master equation coefficients are derived for spectral functions which 
are Laurent series in the frequency with ultraviolet (UV) cut-off A much larger than the 
system frequencies, and infrared (IR) cut-off A much smaller than the system frequencies. 
These coefficients are calculated in detail for arbitrary temperature and arbitrary time after 
the switch-on time. The master equation is solved explicitly for this class of spectra, and 
the general solution for arbitrary coefficients is reduced to solving the classical equation of 
motion of a parametric oscillator. 

It is clearly shown how each coefficient enters into the solution and how the state evolves 
in time. The anomalous diffusion coefficient is actually an "anti-diffusion" term that keeps 
the position uncertainty finite even when the UV cut-off tends to infinity (in contrast with 
the momentum uncertainty). Having the solution gives us a plethora of information such 
as the late time thermal covariance [Eq. (IV. 21)] and uncertainty function (see Sec. IV B 2). 
These results generalize the work of Anastopoulos and Halliwell [15], who already found 
the late time state to be a Gaussian and placed lower bounds on the uncertainty function 
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which we discuss in Sec. IV B 2. This also crystalhzes the work of Hu and Zhang [13] on the 
generahzed uncertainty function for Gaussian states. 

Finally, we derive the master equation that includes the influence of a classical external 
force and solve the master equation in that well. 

B. Systematic Overview 

The paper is organized as follows. We begin our derivation of the master equation 
in Sec. II with a quantum oscillator linearly coupled to a thermal reservoir of quantum 
oscillators. The spectral density here is assumed to be ohmic with high-frequency cut-off A. 
The particle begins both decoupled from and uncorrelated with the reservoir. The coupling 
is turned off at the initial time t = and is switched on within a short time scale r. At 
moderate times, t ^ t (but not necessarily t > fi~^,7^^), we have obtained the master 
equation coefficients as an expansion from zero temperature (Sec. II Bl), as an expansion 
from high temperature (Sec. II B 2), and in closed form as an approximation that is exact at 
both zero and extreme temperatures (Sec. II B 3). Exact closed form solutions are possible in 
terms of various special functions, but they do not straightforwardly reveal the behavior seen 
in our expansions and approximations. Most importantly, we have closed form solutions for 
the master equation coefficients at late times, t ^ 7,^^, and arbitrary temperature (Sec. II B). 

In Sec. Ill we extend these results first to analytic spectral functions (Sec. IIIB) and 
then to Laurent-series spectra (Sec. IIIC). With the inclusion of subohmic terms we must 
introduce a low frequency cut-off A. We take both the UV and IR cut-offs to be very 
large and very small respectively and we only consider spectra that give contributions to the 
master equation coefficients no more divergent than in the ohmic case (which has logarithmic 
dependence on the UV cut-off). This limits our study to equations of motion for the system 
trajectories with the same form as in the ohmic case and whose contributions to the master 
equation coefficients are not markedly different either. Spectral densities of arbitrary powers 
have been studied numerically in the past, but we do not consider them as they correspond 
to classical paths that involve fractional calculus {i.e., integro-differential equations) and 
would be difficult to solve analytically. 

In Sec. IV we solve the master equation for this class of spectra. It is seen that the 
initial solution undergoes damped oscillations while evolving into a Gaussian state of thermal 
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equilibrium. All the cumulants of the Wigner distribution are easily determined as functions 
of time (Sec. IV B 1). In particular we can provide the exact late-time uncertainty function 
(Sec. IV B 2). In Sec. IV C we extend this solution to more general spectra, provided that 
one knows the form of the master equation coefficients. The solution to this master equation 
is reduced to solving for the motion of a classical parametric oscillator. It is seen that simple 
damped oscillations may now be parametric damped oscillations, otherwise the form of the 
solution is the same. 

Finally in Sec. V we extend the method of Calzetta, Roura, and Verdaguer [IS] to include 
the influence of a classical force acting upon the oscillator. The master equation turns out 
to be what one would naively guess. We solve this master equation in Sec. VB and it is seen 
that the external force drives the mean around while the evolution of the variance remains 
unchanged and determined by the reservoir. 

In the last section we conclude with a list of our findings and suggestions for their appli- 
cations. 

II. MASTER EQUATION COEFFICIENTS FOR AN OHMIC SPECTRUM 
A. Ohmic Spectrum Master Equations 

The Lagrangian of a system consisting of a quantum Brownian oscillator with mass M, 
natural frequency Q and coordinate x coupled with coupling constants c„ to an environ- 
ment at temperature T made up of n oscillators with mass m^, natural frequency Un and 
coordinates Xn is given by 

L=^M _ n^x^) + J2 {xl - ivlxl) - d^{t) CnXXn. (II.l) 

n n 

where 9r{t) ^ e"'^*) is a switch-on function with a very short characteristic time-scale r. 
In this section we work with an ohmic bath so that the spectral density function is given by 

J(cu) = -7oMcj (0 < < A), (II.2) 

TT 

where A is a high-frequency cut-off of the bath and 70 is the dissipation frequency. We will 
be using this hard cut-off for all our calculations. Gaussian and exponential decay cut-off 
functions, of the form e"'-'^/^^^ and e"'^/^ respectively, are also common choices. 
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The system and environment begin decoupled at t = 0, but are fully engaged in a time 
r which is assumed to be very short but non-zero. This time-dependent coupling is easily 
dealt with following the work of Hu and Matacz[l 1] on the time-dependent problems where 
all parameters of the system and bath oscillators and their couplings are allowed to be time- 
dependent. When only the system-environment couplings are time-dependent, as in our 
case, the dissipation and noise kernels are given respectively by 

fi{t, s) = 27oM^5A(t - s) er{t)er{s), (II.3) 

u{t,s) = dujcos[uj{t- s)\coi\i(^^^I{uj) er{t)er{s), (II.4) 

where 

e . N 1/"^ , . N sin(At) 

Sj^{t) = - dujcos{tot) = — ^ — (II.5) 

TT Jq Tit 

and it tends to a Dirac delta distribution in the limit A ^ oo. Had we chosen a Gaussian 
high-frequency cut-off function instead of a hard cut-off, we would have found the dissipation 
kernel (II. 3) to be the derivative of a Gaussian distribution that corresponds to an alternative 
representation of the delta distribution in the infinite cut-off limit. On the other hand, the 
integral for the noise kernel (II. 4) is more complicated, even without the cut-off. Although 
it can also be expressed in terms of distributions, for the purposes of this paper it is best 
left to the last. 



1. Initial Time Divergences 

In HPZ[12] the system-environment coupling was taken to be constant in time, yet the 
system and environment were initially uncorrelated. This gives rise to initial "jolts" in the 
normal diffusion coefficient of the master equation with a characteristic time-scale of order 
A^^ and an amplitude proportional to A, which diverges in the limit A oo. It also gives 
rise to a delta-like term at the initial time in the equation of motion: 

u{t) + 27ou(t) + n^u{t) = -47o5A(t)M(0). (II.6) 
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In the limit of infinite cut-off this sudden frequency change gives an initial kick to the 
homogeneous solutions of the equation of motion, which can be seen as follows: 



u = v Au = \im dtu = {11.7) 







i) = -2-fo{t)v - ni{t)u Av = \iml dtv = -2-fou{0) (11.8) 







nl{t) = nl + 4^om (11.9) 

Thus, the classical paths experience a finite velocity change within an infinitesimal time. 
Following the approach in Ref. [ I "^j one can easily see that this kick translates into a distortion 
of the Wigner distribution from the bare initial state to a shifted one 

M^bare(x,p) ^ W.enix , p) = Wi^^^^ , p + 2 M ^qx) . (11.10) 

One expects that this result, which corresponds to replacing 6\{t) with a Dirac delta in 
Eq. (11.6), is a good approximation whenever all the relevant time-scales (f^^^ 70 and t) are 
much larger than A^^. 

The physical origin of the jolts in the coefficients of the master equation as well as other 
initial time divergences (such as the divergent contributions to correlation functions of system 
observables that are due to divergent boundary terms at the initial time (see Appendix D in 
Ref. [17]) can be understood as follows. In general when a system couples to an environment 
with an infinite number of modes, well-behaved states exhibit correlations with arbitrarily 
high-frequency modes. In contrast, states that are uncorrelated for sufficiently high frequen- 
cies (and hence completely factorizable states in particular) are pathological. For instance, 
they have infinite energy and their Hilbert space is even unitarily inequivalent to the space of 
physical states, spanned by the basis of energy eigenvectors of the whole system Hamiltonian 
including the system-environment interaction. (Of course for a finite UV cut-off there are no 
divergences, but the potentially divergent terms are very sensitive to changes in the value of 
the cut-off.) Physically acceptable initial states that correspond to the thermal equilibrium 
state for the whole system can be obtained using Euclidean path integrals [ ]. However, 
the instantaneous preparation functions employed in Ref. [■>] to produce other states in ad- 
dition to the thermal equilibrium state still give rise to initial divergences, as explained in 
Ref. [()]. In order to obtain finite results, one needs to prepare the new initial state within 
a non- vanishing time [7], which corresponds to a physically more realistic situation. The 
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alternative approach that we follow here is to switch on the interaction smoothly within a 
time r much longer than but shorter than any other relevant time-scale of the prob- 
lem. In this way the factorized initial state, which is perfectly acceptable in the uncoupled 
case, becomes adequately correlated with the arbitrarily high-frequency modes in a regular 
fashion. 

When adding the short time switch-on function (11.1) to the spectral density to turn on 
the interaction gradually, the initial jolt is no longer present in the results for the master 
equation coefficients. Moreover, the equation of motion exhibits a smooth transition between 
the decoupled and the coupled system without the term proportional to S/^{t) on the right- 
hand side: 

u{t) + n^u{t) = 0, t < r (11.11) 
u{t) + 2-fou{t) + nlu{t) = 0, T<^t. (11.12) 

However, since the constant bare frequency is of order A so as to cancel out the divergent 
contribution that arises when integrating out the environment and give a finite value of the 
renormalized frequency at late times (larger than the switch-on time), the frequency will 
change significantly in a short period of time of order r. As long as r is much shorter than 
all the other relevant time-scales in the problem (except for A~^), the dynamics for t ^ t 
can be obtained by approximating the time-dependent frequency by a delta function. More 
specifically, since the renormalized frequency changes within a characteristic time-scale of 
order r and with an amplitude of order A, it can be written as 

nl{t) = nl-2cAT5rit), (11.13) 

where fl"^ is the asymptotic constant value for times much larger than r, c is a constant of 
order one, and Sr{t) is a function peaked at t = with amplitude of order and width of 
order r that becomes a Dirac delta in the limit r ^ 0. Therefore, for t ^ r the dynamics 
is governed by 

u{t) + 2-fou{t) + nlu{t) = 2c At S{t)u{0). (11.14) 

Note that At ^ 1 if the switch-on time is sufficiently large to cure the initial divergences 
discussed above. 

The term on the right-hand side of Eq. (11.14) has the same form as in Eq. (11.6). It will, 
therefore, give the same kind of initial kick to the solutions of the equation and generate the 
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same kind of transformation of the reduced Wigner function: 

W^bare(a;,p) ^ Wrenix,p) = Wbare (x, p - C AtMx) . (11.15) 

This phase-space transformation has a Jacobian with determinant equal to unity: 



1 

-cArM 1 



detL = l. (11.16) 



Therefore, it is simple to calculate renormalized expectation values in terms of bare expec- 
tation values and vice versa: 



(A(x,p)),en = JJdxdpA{x,p)W,^n{x,p), (11.17) 
(A(x,p))rcn = (A(x,p + cArMx))bare. (H-IS) 

We can immediately see that the normalization, linear entropy (see Sec. IV B 3) and state 
overlap are all unchanged by the kick. We can also check that the Heisenberg uncertainty 
relation is also preserved as follows. First, we start with the covariance matrix for x and p 
corresponding to the Wigner distribution 



(^xx ^xp 



^px ^pp 



(11.19) 



with a^x = {xx)^cn, o'xp = o'px = {xp)rcn and (7pp = (Tpp = {pp)ren, and which transforms in 
the following way under linear phase space transformations: 

a L^aL . (11.20) 

Hence, from Eq. (11.16) we have 

det (Tbare = det (Tren ■ (11.21) 

Finally, one takes into account that 

(deto-)>^, (11.22) 

corresponds to the formulation in terms of the Wigner function of the generalized Heisenberg 
uncertainty relation due to Schrodinger [8, 9]: 

(AxfiApY - {{x - {x),p~{p)})' > ^, (n.23) 
10 



where {A, B} = AB + BA. In addition to being real and normalized, a phase space distri- 
bution corresponds to the Wigner function of a physical state only when the corresponding 
density matrix (undoing the Wigner transformation) satisfies the condition < p. For 
Gaussian distributions this is guaranteed if and only if Eq. (11.22) is satisfied [1 0]. Thus, 
the kick does not change whether a given Gaussian distribution corresponds to a physical 
state. On the other hand, for non-Gaussian distributions Eq. (11.22) is no longer a sufficient 
condition for that, but one can use the following argument. For a finite switch-on time, the 
frequency change simply corresponds to unitary evolution associated with a time-dependent 
Hamiltonian, which leaves the eigenvalues of the density matrix non-negative. By continuity, 
the limit of a very short switch-on time cannot make these eigenvalues negative. 

If one is interested in studying the evolution of a certain state of the system properly 
correlated with the environment (such as the global equilibrium states considered in Ref. [3] 
or states prepared from those in a finite time) one can always consider the reduced Wigner 
function associated with that state and invert Eq. (11.22) to obtain the corresponding initial 
Wigner function before the interaction was switched on. (Note that for times larger than 
the switch-on time the result is then essentially equivalent to having introduced a time- 
dependent counterterm for the frequency so that the renormalized frequency was constant 
in time.) 



2. Expressions for the Master Equation Coefficients 



The HPZ master equations for the reduced density matrix pr and the reduced Wigner 
function are given respectively by 



d 

d_ 
di 



Hr,Pr + r [X, {p,Pr}] + iD^p ([£, [p, Pr]] + [p, [x, Pr]]) - iDpp [x, [x , Pr]](Jl .24) 

(11.25) 



Wr = {Hr,Wr} + 2T—pWr-2D,p-—Wr + Dpp — Wr, 

op ox op op^ 



where the detailed form of the dissipation function r(t) and the diffusion functions Dxp{t) 
and Dpp{t) can be found in [ ]. Following the derivation of Calzetta, Roura, and Verdaguer 
[18] the coefficients in the HPZ master equation for an ohmic bath with a high cut-off A and 
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evaluated at times larger than the switch-on time are given by^: 

^2 - (11.26) 

(11.27) 



r = 70, 

1 /•* 

Dxp = ds u{t,s)Gret{t,s), 

^ Jo 

/■* d 

Dpp = M J ds u{t,s)—G,et{t,s), 

Grctit,s) = ^sinn(t-s)e-^«(*-^), 



11.28) 
11.29) 
11.30) 



where = a/^^^ " 7o ^^"^ G^etit^s) is the retarded Green function associated with the 
differential equation that results from multiplying Eq. (11.12) by M. According to the 
discussion in the previous subsection, when considering times much larger than the switch- 
on time, its effect on the equation of motion can be approximated by a delta function, as 
seen in Eq. (11.14). Since Gretit-iS) for t < 0, the contribution from the delta term on the 
right-hand side of Eq. (11-14) vanishes when solving for the retarded propagator. Therefore, 
it is indeed independent of the initial kick and given by Eq. (11.14). There is still some 
dependence on the switch-on function 9t-{s) in the noise kernel for s < r, but for t^r the 
contribution from s < r to the time integrals in Eqs. (11.28)-(11.29) is negligible provided 
that the behavior during the switch-on time is sufficiently regular, which is indeed the case. 
Similarly, the behavior of the master equation coefficients during switch-on time, which we 
do not calculate here, will be sufficiently regular so that its contribution to the evolution of 
the reduced Wigner function (or density matrix) is also negligible when considering times 
larger than r and one only needs to include the effect of the kick on the initial Wigner 
function, given by Eq. (11.15). 

Some clarification is in order here about the use of a finite frequency cut-off A for the 
environment spectrum. Having a finite cut-off in the frequency integral for the noise kernel, 
given by Eq. (11.4), is necessary to obtain a finite result for the coefficients of the master 
equation because, as it will be seen below, it gives rise to contributions to the diffusion 
coefficients that become logarithmically divergent in the limit of large A. In contrast, the 
use of a finite cut-off in the frequency integral of the dissipation kernel is not required to 



^ This method has the nice property that the master equation coefficients after the switch-on time are 
exphcitly independent of any reasonable behavior during the initial switch-on. 
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obtain a finite result. In fact, for computational convenience we are going to take A ^ oo in 
this case because then the dissipation kernel becomes local and we just have to deal with an 
ordinary differential equation rather than an integro-differential one. One may object that 
treating the noise and dissipation kernels on a different footing could lead to inconsistencies 
{e.g. the fluctuation-dissipation relation will not hold exactly any more, as pointed out by 
HPZ [ ] with regard to the treatment of this problem in Ref. [■>]). However, our results are 
consistent when properly understood. If one considered an expansion of the exact master 
equation coefficients in terms A, one would have leading terms of order log A, terms of order 
one and subdominant terms involving inverse powers of A. Our approach gives the right 
result for the log A and order-one terms: using the local approximation for the dissipation 
kernel only alters the result for the subdominant terms. This can be seen as follows. When 
integrating 6\(t — s), given by Eq. (II. 5), with a sufficiently well-behaved function /(s), one 
has 

/oo 
ds5A{t - s)f{s) = fit) + f'it)A-' + 0(A-2). (11.31) 
-oo 

This condition will indeed be satisfied in our case since all the other relevant time-scales in 
the problem are much larger than A~^ and all the functions that 6\{t — s) will be integrated 
with are regular and extremely uniform over that scale. Therefore, the discrepancy between 
using the Dirac delta function and using 6^(1 — s) corresponds to terms of order A~^ or 
higher inverse powers, which give contributions to the master equation coefficients of order 
A^^ log A or higher. 

In order to compute the coefficients of the master equation, we first perform the relatively 
simple time integrals in Eqs. (II.28)-(II.29). On the other hand, the frequency integral in 
the noise kernel, as given by Eq. (II.4), is particularly difficult due to the 1/ sinh(ci;/2T) in 
coth(tt;/2T). It is this integral that we save for last. Others[ i, ] have also seen the utility 
or at least the simplicity of reducing the master equation to a collection of one-dimensional 
frequency integrals. We are able to obtain four sets of new results: (a) exact results at late 
time for arbitrary temperature, (b) asymptotic expansions for high and low temperatures, 
(c) approximate general results for arbitrary temperature and all times (larger than the 
switch-on time), (d) time-dependent diffusion functions in closed form at high temperatures 
for all times (larger than the switch-on time). 

After computing the time integrals in Eqs. (11.28) and (11.29), the diffusion coefficients 
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become 



D,,it) = +^ (FI3 - n^Fh) 

cos(Ot)e-^°* (FC3(t) - n^FC.it) + 27oFS2(i)) 

TT 

+^ sin(!1i)e-'" (i„ (FC3(«) + fijFCiC*)) + ((1' " '!«) FS^M - FS,{t)'j , (11.32) 



TT 



cos(f2t)e-^«* (27oFC3(t) + l^'FS2(t) - FS^{t)) 
sin(f]t)e-^°* (-fi^FCi(t) + (q' - j^) FCs{t) + 70 (nlFS^it) + FS.it)) 



(11.33) 



in terms of a single family of frequency integrals FCat, FSat, FI 



N'- 



UJ 



N 



coth (^) 



FItv = / , duj, (11.34) 







Fc^w = / : : z:r rry u, (11.35) 



cj^' cos(a;t) cotli(^) 

sin(a;t) coth (^) 
uj^ - nif + 472cu2 
FI^ = FC^(O), (11.37) 

FS^(t) = -^FC;v-i(t), (11.38) 
FC^(t) = +^FS^_i(t). (11.39) 

Fortunately the only integral that needs to be computed is FCi(t). All other integrals can be 
generated from this one integral. In general the integrals cannot be solved directly, therefore 
it is necessary to expand the hyperbolic cotangent into a series of simpler functions. This 
will be done with low and high-temperature expansions. 

When computing the time integrals that led to Eqs. (II.32)-(II.33), we just took the 
switch-on functions equal to one, so that sufficiently simple analytical results could be ob- 
tained. This means that, within our approach, the results in Eqs. (II.32)-(II.33) are only 
valid for t ^ r. In fact, for t < r, they essentially coincide with those of Ref. [12] for 
an ohmic environment. For instance, the term FS4(t) in Eq. (11.33) gives rise to the same 
initial jolt, with a width of order and an amplitude proportional to A, found in Ref. [12]. 
However, as we discussed above, the results should be valid for times larger than the switch- 
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on time because the contribution from the switch-on period to both the master equation 
coefficients and the evolution of the density matrix at t ^ r is neghgible. 

It is interesting to note that the coefficients -D^p(t) and Dpp{t) both exhibit logarithmic 
divergences in the limit A oo (for times larger than the switch-on time and thus larger 
than A~^) due to the term proportional to Fl^."^ This has been pointed out for Dxp{t) in 
Ref. [11], where the coefficients of the master equation were calculated perturbatively to 
second order in the system-environment coupling constants (linear order in 70). The fact 
that there is also a logarithmic divergence in Dpp{t) was not seen in that reference because 
it is quartic in the system-environment coupling constants (quadratic in 70), as it can be 
seen in Eq. (11.33). Moreover, such kinds of perturbative calculations cannot be employed 
to study the long time behavior since they are only valid for 7ot -C 1 and they miss for 
instance the exponential decay of the second and third terms on the right-hand side of 



B. Exact Late-Time Behavior and Approximate High and Low- Temperature Re- 
sults 

For the late-time diffusion coefficients, only two integrals need to be performed and they 
can be expressed in closed form with reasonably intuitive functions:^ 



where terms involving negative powers of the cut-off have been neglected and H(z) is the 
harmonic number function defined in Appendix A 1. It behaves like log (z) except near the 

^ As it is well-known and can be checked from Eq. (11.34), the term logarithmic in A is not present in FI3 
when considering a fixed finite cut-off A and temperatures much higher than A. 

^ Many of the expressions derived throughout this paper assume underdamping with an extreme cut-off, 
i.e., 7o < < A with fl = -^/fi^ — --y^. They can be used for the overdamping regime by making the 
following analytical continuation: f2 — > 17 with 7 = •\/7q — fi^ real. 

Therefore, Eqs. (II.40)-(II.41) can be applied to the overdamping case if the Im and Re terms are first 
expanded assuming Ct is real, and then the analytical continuation ^ 17 is made. 



Eqs. (II.32)-(II.33). 




(11.40) 



(11.41) 
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origin, where it does not diverge and actually vanishes. The Re[- ■ ■ ] term in FI3 is effectively 
the \og{A/Qr) divergence. It disappears at extreme temperature (higher than the cut-off). 



1. Expansion from Zero Temperature 



A low-temperature expansion of the master equation coefficients can be obtained using the 
following expansion for the hyperbolic cotangent in the corresponding frequency integrals: 



coth {^) 



1+2J2 



(11.42) 



k=l 



Moreover, since the results for all the relevant integrals can be obtained from FCi(t) simply 
by differentiating with respect to t, we will concentrate on computing Eq. (11.35) with = 1. 
In order to do that, one first decomposes the integrand in Eq. (11.35) into partial fractions 
and splits the integral duj into duj — dw. Finally, using Eq. (11.42) they become a 
linear combination of exponential integrals. The result can be written as 

1 El f+(7o - ih)t] El f-(7o - iVL)?^ 



FCi(t) 



47ol] 



— Im 



-(70-«n)t 



+ 



TT 



^ cos(nt)e-^«* + AFCr(t) + AFCi^^(t) 



3+(70-«f^)i 



470^^ 



(11.43) 



where Ei(2;) is the exponential integral detailed in Appendix A 2. The first three terms 
come from the f^^du integral of the constant term in Eq. (11.42). On the other hand, 
AFCf"'"(t), which vanishes at zero temperature, is the contribution from the dtu integral 
of the remaining terms in Eq. (11.42): 

Finally, the term AFC^'^(t) contains all the dependence on the cut-off A, which arises from 
all the doj integrals. In particular, its value for zero temperature, which corresponds to 
taking only the first term in Eq. (11.42), can be calculated explicitly: 



AFCf''(t) 



870^^ 



— Im 



El ( +(70 - i{n + A))t j El ( -(70 - + A))t 
El (+(7o-z((]-A))t 



+ 



El 



-(70-«n)t 



g+(70-«n)t 
g+(70-«f^)i 



(.11.45) 
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Note that, since FCi(t) is finite in the hmit A oo, the contribution from AFC^^(t) will be 
very small for sufficiently large A. However, when differentiating several times with respect 
to t, this term gives rise to the divergences that FC3(t) and FS4(t) exhibit during a short 
initial period of order A~^, which are responsible for the initial jolt of the master equation 
coefficients found in Ref. [12]. Moreover, even in those cases it would be enough to consider 
just AFCi'^(t) rather than AFC\^{t) because the integral in Eq. (11.44) is convergent for 
arbitrary positive powers of the frequency due to the exponential factor. 

In fact, given that the divergent behavior for short times is irrelevant when considering a 
switch-on time sufficiently longer than A~^, it turns out that in our case we can completely 
neglect the contribution from the cut-ofT dependent terms for sufficiently large A. It would 
only be necessary when obtaining FI3 by evaluating FC3(t) at t = 0, but it is much simpler 
and more accurate in general to use the exact result for FIi and FI3 computed earlier and 
given by Eqs. (II.40)-(II.41). 

To sum up, taking into account the remarks in the previous paragraph together with the 
fact that the contribution from the third term in Eq. (11.43) cancels out when adding all the 
time- dependent terms in Eqs. (II.32)-(II.33), all that one needs [in addition to Eqs. (11.40)- 
(11.40)] to calculate the master equation coefficients is 

El (+(70 - zn)t) El (-(70 - 2fi)t' 



AFCi(t) = -^Im 



g-(70-«f1)i ~'~ g+(7o-«n)t 



/ ^^ :—du}. (11.46) 

Given the asymptotic behavior of Ei(z) for large l^l, one can see that the Im[- ■ ■] terms 
decay only weakly in time like l/{VL^tY\ they both have complicated VL frequency oscillations. 
The integral in the last term, which vanishes when the temperature tends to zero, can be 
performed explicitly, but the sum becomes particularly complicated. 



2. Asymptotic Expansion from High Temperature 

A high-temperature expansion can be obtained by making use of the following expression 
for the hyperbolic cotangent: 

/ \ 2T ^ ATuj 
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which will yield a series solution that is best in the high temperature, late time regime 
27iTt ^ 0. FCi(t) can then be calculated in several steps. 

First, one calculates the du integral corresponding to Eq. (11.35) considering only the 
first term in Eq. (11.47), decomposes the resulting integrand into partial fractions and ends 
up with a linear combination of Exponential Integrals. After simplification, the result can 
be written as 

FCf^(t) = 7^ (n cos(fit) + 7o sm{nt) e"^"*. (11.48) 

Next, we perform the du integral corresponding to Eq. (11.35) with the general term in 

the sum of Eq. (11.47) replacing the hyperbolic cotangent. After decomposing into partial 

fractions we end up again with a linear combination of Exponential Integrals, and the final 

result after simplification is 

°° 7^-70* r 
FCf (t) = y y — ^ -4A;7rTf27oe(^°-2'='^^)* 

+ 2 {AkVT^ + 7^) fi2 + (^2 _ 4k2^2rp2f\ L 

+70 (^-Ak'^n'^T'^ + + 7o) sin(fit) + (^Ak\^T^ + + 7o) cos(fit) . (11.49) 

The sum over k can be performed explicitly for the terms multiplying the sine and the cosine, 
so that we get 



^ sinh p cos(fit) + sin(^)sin(fit) 

FCf (t) = e--* + AFEf (t), (11.50) 

470^^ cosh f p j - cos 

with 

AFELT(t) = ' y (11.51) 

which satisfies the following property: 

^AFE^f (t) = -(27rr)AFEjf^,(t). (11.52) 

Putting everything together, we have 

^ sinh (^) cos(nt) + sin (^) sin(nt) 

FCi(t) = ^ — e-^«* 

470^ cosh (^p) - cos (f) 

+AFELT(t) + AFCf^(t), (11.53) 

where AFCf^(t) contains all the dependence on the cut-off A, and comes from all the duo 
integrals. Nevertheless, for the same reasons given in the previous subsection, in our case one 
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can neglect this term for sufficiently large values of A. Furthermore, since the contribution 
from the first term in Eq. (11.53) cancels out when adding all the time-dependent terms in 
Eqs. (II.32)-(II.33), and provided that we use Eqs. (II.40)-(II.41) to compute FIi and FI3, 
it is sufficient to consider 

1 °° h.p-2-KTtk 

AYCAt) = -, — y (11.54) 

(-2vrT)^^(p+(n^)^)_4(^)^p' 

in order to calculate the master equation coefficients. 

All terms on the right-hand side of Eq. (11.54) decay with temperature-dependent time- 
scales. Since it cannot be expressed in closed form with any intuitive functions, the sum 
is not explicitly performed. Although useful for numerical evaluation, this form does not 
immediately reveal the 7ot and f2t behavior. Therefore, it is convenient to approximate the 
sum by an integral that can be computed explicitly, as follows: 



fe=i -^^ 

Ei (^(70 27rT)t^ Ei (^(-70 + tQ - 27rT)t j 



AFCi(t) ^ ^^Im 



. (11.55) 



From this expression one can see that the most lingering terms decay exponentially as 
e"^'^"^* (9(l/t). There are also oscillations with frequency Cl. 



3. Approximate General Solution 



By inspecting the corresponding terms in the high (11.55) and low-temperature (11.46) 
expansions, one can get a good idea of the general solution: 



FCi(t) 



TT 



sinh ( p ) cos(fit) + sin (^) sm{nt) 



-lot 



470^^ 
1 



cosh ( p 



— Im 



47ol] 
-AFCf^(t). 



El 



cos(^) 

7o + «fi + 2TiT)t I El ( (70 + 27TT)t 



(11.56) 



This expression is exact at both zero temperature and high temperature, and it has only mi- 
nor discrepancies within a small range of intermediate temperatures. Moreover, as explained 
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above, all that one really needs to calculate the master equation coefficients is 
1 



AFCi(t) 



— Im 



El ( (-70 + zfi + 2TTT)t ) El ( (70 + 2TiT)t 



(11.57) 



4- Comparison with Numerical Results 

From our analysis of the integrals it can be seen that there is only one fluctuation fre- 
quency Vt and it is always damped. There are only two damping rates, 70 and 27rT. All 
temperature-damped terms also have coupling damping, so one can say that t 70"^ is 
always late time, although the temperature can hasten this. 

Since we have solved the integrals in the different temperature regimes, we will now plot 
the lowest order terms of each expansion in the different temperature regimes. The coupling 
will be set at a moderate level to enhance their differences. 




Qrt 



FIG. 1: The case of a moderate temperature. FCi calculated — numerically, to first order in the 
• • • low-temperature expansion, to first order in the • • • high-temperature expansion, and with 
the — approximate general solution. All of the approximate solutions are fairly close with some 
minor discrepancy at short times. The low-temperature expansion will continue to be slightly off 
into late times. 



Note that the high-temperature expansion (including the first order correction) gives 
very good agreement even at much lower temperatures T ~ fi,. (much lower than the regime 
T > A where the high-temperature expansion is clearly expected to be very accurate). 
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FIG. 2: The case of a high temperature. The — approximate general solution and the • • • first 
order high-temperature solution stick very well to the — numeric solution, while the • • • 
first order low-temperature solution has amplitude and phase discrepancies. 




FIG. 3: The case of a low temperature. The — approximate general solution and the • • • first order 
low temperature solution stick very well to the — numeric solution, while the • • • first order 
high temperature solution has some amplitude and phase discrepancy initially. 

When we say n^^ order in the high-temperature expansion, we more acurately mean up to 
the e"^'^-^*"' term in the high-temperature expansion. Strictly speaking this is not a temper- 
ature expansion but an expansion in e"^'^"^*, the only terms which decay at a temperature- 
dependent rate, so it is really a high-temperature, late-time expansion. Similarly, when we 
say n^^ order in the low-temperature expansion, we mean keeping the n^^ order term in the 
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sum appearing in the last term of Eq. (11.46). 



C. Some Closed Form Results 

1. Extreme Temperature, Arbitrary Time 

Considering just the first term on the right-hand side of Eq. (11.47) for the hyperbohc 
cotangent, one gets the following result for temperatures much higher than the cut-off: 

T > A, 

D^p = 0, (11.58) 
= +2M7or. (11.59) 

This is the common, but questionable high-temperature limit. There is no first order term 
for the anomalous diffusion coefficient, but there are lower order terms that vanish only for 
infinite temperature. They include the term that gives rise to the log A divergence when the 
temperature becomes much smaller than A. 



2. Arbitrary Temperature, Late Time 



Neglecting the terms multiplied by the factor e '^^^ in Eqs. (II.32)-(II.33), we obtain the 

following result, which is valid for times much larger than the relaxation time {i.e., t ^ 7(7^): 

1 



t > 



D 



xp 



7o 

-llm 

Tin 



+ ^Re 

TT 



H 



7o + 
2nT 



H 



7o + 
2nT 



D. 



pp 



+2M7„r+^(f!^-7„2)Im 



TT 



H 



H 



7o + 
2nT 



A 



H 



7o + 



(11.60) 



(11.61) 



2nT J \ 2nT 

The Re[- ■ ■] terms become log{A/Qr) at zero temperature and vanish in the infinite tem- 
perature limit. On the other hand, the value of the Im[- ■ ■] terms is cos^^(7o/^^r) at zero 
temperature, whereas it also vanishes in the infinite-temperature limit. 
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D. Comparison with Caldeira et al. 



We now compare our results to CCR's weak coupling master equation [4] , which differs 
from the HPZ equation and has the following diffusion coefficients: 



D 



c 

xp 



0, 



Dg = 7oM(],coth(^V 



(11.62) 
(11.63) 



The CCR master equation is frequently used for its simplicity and is believed to be accurate 
at late times for weak coupling and arbitrary temperature. One can see from Fig. 4 that, 



M Qj,^ 



100 




0.2 0.4 0.6 



T 

Qr 



FIG. 4: Late time Dpp for — high temperature, ■ • • Caldeira, — HPZ at A = 10^ ftr, and • • • 
HPZ at A = 10"^^. 

ignoring the contribution from the cut-off, the CCR approximation matches extremely well 
with our exact results for the normal diffusion coefficient at weak coupling. The CCR error 
from neglecting the cut-off dependence is determined by the order of magnitude of the cut- 
off scale since it is logarithmic in A. In any case, the CCR approximation underestimates 
the magnitude of the diffusion coefficient. 

Unlike the normal diffusion coefficient, the difference is much more severe for the anoma- 
lous diffusion coefficient, which is completely absent in CCR. The largest contribution (in 
the weak coupling regime) to the anomalous diffusion coefficient comes from the cut-off and 
it does not vanish at finite temperature (see Fig. 5). Moreover, it cannot be regarded as 
vanishing at weak coupling because it is only proportional to one power of the coupling 
constant, which is the order to which CCR's master equation should be valid. 
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FIG. 5: Late time Dxp for — high temperature or equivalently Caldeira, — HPZ at A = lO^flr, 
and • • • HPZ at A = 10^ fir- 

The final point that should be made is that OCR's arbitrary temperature master equation 
is only valid in the weak coupling regime for late times and, unfortunately, late time always 
means t ^ 7q"^ (according to our exact results for the HPZ master equation). Therefore, 
the weaker the coupling, the longer one must wait for OCR's master equation to be valid. 

III. MASTER EQUATION COEFFICIENTS WITH LAURENT SERIES SPEC- 
TRA 

A. Derivation of The General Integrals 

We start by considering the independent effect of power-law terms in the spectral function 
with ultraviolet cut-off A and infrared cut-off A whenever necessary. We take these cut-offs 
to be extremely large and extremely small respectively, not because that should necessarily 
be the case for any relevant physical situation, but to restrict the model to something that 
we can solve explicitly. Hence, we take a spectral function of the form 



where N is an integer number and ooc is a product of characteristic frequencies chosen 
so that well-behaved results are obtained for the master equation coefficients in the limits 
A ^ oo and A ^ 0. Traditionally Uc has been chosen to be A whenever a UV cut-off was 




(A < < A), 



(IH.l) 
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necessary; we will end up making a similar choice here. 



1. The Classical Trajectories 

Given a general spectral function the equation of motion for the system trajectories 
corresponds to an integro-differential equation: 

= u{s) + n\{s) + ^ [ ^{s - s')u{s')ds' . (III.2) 

It is thus convenient to perform a Laplace transform 

/•CO 

/(C) = / e~<f{s)ds, (III.3) 







under which the equation becomes purely algebraic, since the integral involving the dissipa- 
tion kernel is merely a convolution. We then have 



or, equivalently. 



where uq, uq are the initial conditions for u{s). Given the expression of the dissipation kernel 

fi{s) = — / In{u) sm{ujs)duj, (III. 6) 

one can easily compute its Laplace transform: 

o /"A 2 / \ 



m = -M7^_i / — - - du. (III.7) 

71" Jx UJ^ + C^ \UJcJ 

The results for different integer values of are listed in Table I. The column labeled with 
/ic(C) corresponds to the outcome of the integral in Eq. (III. 7) after taking the limits A ^ cxo 
and A ^ wherever that gave a finite result. The values of ujc were chosen so that the 
integrals that appear when calculating the diffusion coefficients were finite (more details can 
be found in the next two subsections). Finally, the column labeled with fidC) follows from 
substituting uoc in fidO with the values of the previous column and taking the limits A — cxd 
and A ^ whenever they give a finite result. 

Had we chosen fractional powers of the frequency in the spectrum or even a different 
functional dependence such as a logarithm, this would have resulted in an integro-differential 
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N- 1 


/ic(C) 




A(C) 


—n 









-4 


(-i^i + ff) 


X^UJc 





-3 


fM7-3(-^logf) 







-2 


M7-2 (-f) 


XuJc 





-1 


fM7_i (-c^elogf) 


UJc 


fM7_i (-c^elogf) 



1 
2 
3 


M70 (C - f A) 
fM7i (?logf-|^) 

fM73(-Slogf + ||C^-i§) 


1 

A 

A3 


M70 (C - f A) 
-fAf7if 
-fM72f 


n 




A" 





TABLE I: Laplace transform of the dissipation kernels for various power-law spectral functions 
after discarding 0{X) and 0{1/A) terms. N is the power of the spectrum: = 1 is ohmic, > 1 
is supraohmic and A^ < 1 is subohmic. ficiC) is the Laplace transform of the dissipation kernel 
before a reasonable choice of characteristic frequencies coH^ ^' is made, whereas is the choosen 
physical dissipation kernel. Constant terms are renormalizations to the bare frequency. 

equation for the classical system trajectories. In addition to being less classical, they would 
also be much more difficult to solve. 

Integer subohmic spectra, (corresponding to — 1 < in Table I) all exhibit a similar 
behavior, except for the = case. For A^ < it is never possible to change the classical 
ordinary differential equation (ODE) while keeping the diffusion coefficients finite in the 
limit A — s> 0. Indeed, the only way to guarantee that all the frequency integrals involving 
the noise kernel that appear when calculating the diffusion coefficients are finite is by taking 
uj^ ~ A"~^ (higher powers of A are allowed, but would give vanishing results for the diffusion 
coefficients). In that case, there are no additional contributions to the equation of motion 
left. In contrast, for A^ = the diffusion coefficients only exhibit a logarithmic IR divergence 
(proportional to log A) when considering a finite non- vanishing value of Uc independent of 
A. Since we tolerated logarithmic dependences on the UV cut-off A in the ohmic case, it 
would be natural to allow a similar situation with the IR cut-off A. However, in that case 
the Laplace transform of the dissipation kernel exhibits a logarithmic term that makes the 
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equation of motion for the system trajectories (which corresponds to an integro-differential 
equation) much more difficuh to solve. Therefore, we will not consider the = case, for 
the same reason why we did not consider fractional powers of the frequency in the spectral 
function or even a non-logarithmic dependence on the frequency. 

On the other hand, integer supraohmic spectra (corresponding to — 1 > in Table I) 
all exhibit a similar behavior. None of them are capable of changing the nature of the 
classical ODE while keeping the diffusion coefficients finite in the limit A ^ cxd. Indeed, the 
only way to guarantee that all the frequency integrals involving the noise kernel that appear 
when calculating the diffusion coefficients are finite is by taking cj" ~ A*^ (higher powers of 
A are allowed, but would give vanishing results for the diffusion coefficients). In that case, 
the only contributions to the equation of motion that are left are frequency renormalization 
terms. 

B. Analytic Spectra Master Equation Coefficients 

In this subsection we consider spectral functions that are analytic in the frequency (and 
vanish in the limit of zero frequency) . They have the following form: 

/|,l(c.) = ^M^7(^) (c.<A), (III.8) 

where 7 is an analytic function. Eq. (III. 8) can be written as 

/|,}(a;) = -M^5^7n(J (^<A), (III.9) 

n=0 

which corresponds to a linear combination of terms with A^ > 1 among those considered in 
the previous subsection. 

1. The Classical Trajectories 

The classical equation of motion associated with these spectra is an ODE that differs 
trivially from the ohmic case: 

C«o + u, = (2^(0 + M70C - -M ^ -^A u{0 + ^]2«(C), (iii.io) 

\ n=0 ^ / 
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which can be written in the more compact form 



(uo + uo= (C' + 270C + ^l) HO 



(III.ll) 



by introducing the renormahzed frequency 



nl = n^- -A 



n + l 



In 



(III.12) 



We can see that only the power (from the ohmic-hke term) gives rise to the dissipa- 
tion term,^ whereas all the remaining supraohmic terms merely contribute to an additional 
renormalization of the bare frequency. In fact, one could rig supraohmic anticoupling terms 
to eradicate the divergent difference between the bare and renormalized frequency, although 
there is no physical motivation to do so. 

2. Master Equation Coefficients 

For times longer than the switch-on time, when our explicit expressions for the diffusion 
coefficients can be used, the only difference between the master equation for an analytic 
spectrum and the ohmic spectrum lies in a modification of the contribution from the FI3 
integrals to the diffusion coefficients. This can be seen as follows. When considering a term 
proportional to {uj/uJc)"' in the spectral function, instead of the FI3 integral of the ohmic 
case, one gets 



and the only way to get a finite non-vanishing contribution is by taking Uc proportional to 
A (which was the choice already made in Eqs. (III.8)-(III.9)). Therefore, one just needs to 

^ Note that, since an additional term uq is missing, the inverse Laplace transform of 27oC'u(C) gives not only 
the local dissipation term 2"fQu{t), but also a term proportional to S(t)u{0), which corresponds to the term 
on the right-hand side of Eq. (II. 6) and is responsible for the initial kick when no smooth switching-on 
function is present. 




(III.13) 



(III.14) 



(III.15) 
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introduce the following simple substitution in the expressions for the diffusion coefficients 
of the ohmic case: 

Ohmic — > Analytic 

f^^O An ^^70 

where we took cUc — and discarded terms involving negative powers of A. Note that 
all other integrals contributing to the diffusion coefficients will exhibit lower powers of the 
frequency in the numerator so that their result will be proportional to inverse powers of A 
and can be neglected. The integrals that correspond to FC3(t) and FS4(t) in the ohmic case 
also exhibit divergences that would lead to a non-vanishing result (and even a divergent 
one for FS4(t)) when dividing by A", but that is only for short times of order A~^ after 
the initial time. This is the regime where our explicit expressions are not accurate if the 
interaction is gradually switched on with a much longer characteristic time-scale, in which 
case the contribution to the diffusion coefficients from that period is very small. 

We close this subsection by pointing out that, since FI3 is where the log A divergence arises 
for the ohmic case, one could rig an infinite number of finite supraohmic anticoupling terms 
to cancel the divergence, though there is no physical motivation to do so. One could even 
ask whether we can both renormalize the diffusion coefficients and keep the bare frequency 
finite and unmodified, which amounts to requiring the following two conditions to hold 
simultaneously: 

00 

E;^--To. (111.17) 

n=l 

00 . 

-70 log ^. (III.18) 

n=l 

The answer is obviously in the negative, at least not with finite couplings. 



C. Laurent Series Spectra Master Equation Coefficients 

In this subsection we extend the form of the spectral function considered in the previous 
subsection to that of a Laurent series. This is done by adding an analytic function of u!~^ 
to the analytic function of cu already considered there, so that we have 

/{,,^}H = /{,}(a;)--M-(^(-) (A<u;), (III.19) 
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where is an analytic function. Eq. (III. 8) can be written as 



n=0 



•n 




(III.20) 



which corresponds to a hnear combination of terms with < among those considered in 
Sec. Ill A. 

As we saw in Sec. Ill A, in the limit A ^ there is no modification of the equations of 
motion due to the terms in Eq. (III. 20) with negative powers of uj. Hence, the dissipation 
and renormalized frequency terms, r(t), r2^g„(t), in the master equation are exactly the same 
as with ohmic and, more generally, analytic spectra: r(t) = 70, fi^en — ^r- However, there 
will be some non-trivial contribution to the diffusion coefficients that we analyze below. 

For times longer than the switch-on time the only difference between the master equation 
for the Laurent-series spectrum and the analytic spectrum lies in a modification of the 
contribution from the terms most sensitive to IR divergences, namely FIi and FCi(t). This 
can be seen as follows. When considering a term proportional to {uJc/oj)'^ in the spectral 
function, instead of the FCi(t) integral of the ohmic case, one gets 



where we used the expansion coth(co'/2T) = (2T/uj) + 0{uj^) in the first equality. The only 
way to get a finite non-vanishing contribution is by taking uj^ proportional to A (which was 
the choice already made in Eqs. (III.8)-(III.9)). While FIi can be obtained by evaluating 
FCi(t) at t = 0, all the remaining integrals appearing in the expressions for the diffusion 
coefficients will involve a less negative power of the frequency and, when dividing them by 
A", will give a result corresponding to positive powers of A, which can be neglected in the 
hmit A 0. 

Thus, in order to obtain the diffusion coefficients for a Laurent-series spectral function of 
the type considered in this subsection, one just needs to introduce these simple substitutions 




(IH.21) 
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in the expressions for the diffusion coefficients of the ohmic case: 



Ohmic — > Laurent 

7oFl3^7oFl3 + ^ ^ = "E-' ("1-22) 



OO 



n=l 

7oFCi(t) ^ 7oFCi(t) - ^ cos (At)0 = ^ (III.23) 

n=0 

7oFIi ^ 7oFIi - ^0 . 

In summary, the supraohmic terms shift the bare frequency and FI3 integrals, and can be 
used to shift their values in the ohmic theory, whereas the subohmic terms shift the FCi 
and FIi integrals. 



IV. SOLUTIONS TO THE MASTER EQUATION 



A. Solutions to the Laurent-Series Spectra Master Equation 



1. Matrix Representation 



First, we express the master equation 



in a more compact form: 



d_ 
di 



Wr = (V^qD Vq + V^Hq) Wr , 



(IV.l) 



(IV.2) 



with 



-if 



D = 



d_ 
dx 

d_ 
dp 



-D 



xp 



Dxp Dpp 



(IV.3) 



(IV.4) 



This is a hyperbolic second order partial differential equation (PDE). As D is a function of 
time, in general the equation is not seperable in time. It is not seperable in phase space 
either. 
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2. Phase-Space Fourier Transform: The Characteristic Function 



The nature of the PDE suggests a Fourier transform of the phase-space variables since 
derivatives are more comphcated than algebraic parameters. Furthermore, not only does 
a Fourier transform reduce the PDE to first order, but the computation of expectation 
values also becomes trivial since we are then working with the characteristic function of the 
distribution. 

The Fourier transform is defined as 

/oo i*oo 
dx dpe-^'^-V(q), (IV.5) 
-oo J ~oa 



and it exhibits the usual properties: 

The master equation then becomes 

d 



f r\ poo poo 

^(0) = / dx dp g;/(q). (IV.6) 

i J —oo J —oo 



^^Wr = {tk^Btk + zk^mVk) Wr, (IV. 7) 

where Wr = JFjPI/j,} and the normalization of Wr{t,q) implies Wr(t,0) = 1. Finally, it is 
convenient to group all the derivatives on the left-hand side, so that we have 

d 



, k^HVk ) Wr = -k^Dk Wr. (IV.8) 



3. Method of Characteristic Curves 

The method of characteristic curves involves looking for parameterized curves in the 
domain (t, k) along which the first order PDE becomes a set of first order ODEs. For each 
one of those curves we have 

Wr{t,k) = Wr{t{s)Ms)), (IV.9) 

d dt d (ik^ 

= :r7r,^r + -r ^^w. , iv.io 

ds ds at ds 

Next, we attempt to match the right-hand side of Eq. (IV.IO) to the left-hand side of 
Eq. (IV.8). This results in a system of ODEs in the parameter s. We will look for curves 
that synchronize with the initial time so that t(0) = 0, k(0) = kg. The time solution is 
simple: 

^ = 1 ^ t(s) = s. mil) 

ds 
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On the other hand, the solution for the Fourier transform of the phase-space variables is a 
bit more involved: 



ds 



= k^e 



In order to calculate the exponential of the H matrix we diagonalize it, so that 



k^Hq± = h±k^q±. 



The eigenvalues, eigenvectors, and exponential matrix are given by 



(IV.12) 



(IV.13) 



h± = 7o±^^^ q± 



1 

-Mh± 



-tu 



2M70 + M{d/dt) 



-M'^Q.l M{d/dt) 



Gret{t, 0), 



(IV.14) 

(IV.15) 



where GT:et{t,0) is the retarded Green function in Eq. (11.30). We now have the rules for 
transforming back and forth between the domain coordinates (t, k) and the characteristic 
curve coordinates (s, ko). kg does not change along the characteristic curve, but for a given 
s it uniquely specifies a particular curve (except at the origin, where all curves intersect). 

Using these results, we can immediately apply the method of characteristic curves to 
solving Eq. (IV. 8) as follows: 

d 

ds 
d 

ds 



-Wr{t{s),k{s)) = -k^D(t)k W,(t(s),k(s)), 
;>V. (s,e^"^ko) = -k^e^"D(s)e^^ko W.(s,e^^ko 



(IV.16) 
(IV. 17) 



The last equation is a linear ODE whose solution can easily be found to be 
W.(.,e^«"ko) =W.(0,ko)e-/o'^^''^oe='"°(^')e='""''o, 



(IV.18) 



where Wr (0, ko) is the initial reduced Wigner function at t = 0. We can now express the 
solution back in terms of k and get the final result 



Wr (t, k) = Wr (0, e-*«"k) e-^k^'^rWk^ 



with 



<7 T 



{t) = 2 re(^-*)HD(s)e(^-*)^ds. 
Jo 



(IV.19) 



(IV.20) 
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One can check that this result agrees with those in Ref. [ I S] if one considers equal times in 
the two-point correlation functions there. 

From Eqs. (IV. 15) and (11.30) one can see that all the integrals in Eq. (IV. 20) involve a 
factor e^'-*"*-'''''' times some oscillatory factor. Taking into account that the diffusion coeffi- 
cients tend to constant asymptotic values Dxp and Dpp for sufficiently late times, one can 
see that for t ^ 70 the integral is dominated by large values of s. Therefore, the asymp- 
totic value of (TT{t) can be calculated using the constant asymptotic values of the diffusion 
coefficients in Eq. (IV. 20) and taking the limit t 00, with the following result: 



The solution (IV. 19) has two factors. The first one tends to one in the long time limit 
and encodes the disappearance of the initial state (we will call it the death factor). The 
second factor describes the appearance of a Gaussian state that evolves in time and tends 
asymptotically to a state that corresponds to thermal equilibrium (we will refer to this as 
the birth factor). All initial distributions evolve towards this final Gaussian state, whose 
covariance matrix is given by Eq. (IV.21). This state does not look like the thermal state 
of a free harmonic oscillator because of the coupling to the environment. It results from 
considering the thermal equilibrium state for the whole system (system plus environment) 
including the system-environment interaction, which gives rise to a non-trivial correlation 
between them, and tracing out the environment. 

The death factor contains the information on the initial conditions, describes the gradual 
disappearance of the initial distribution and it is always temperature independent. The 
initial distribution undergoes damped oscillations with characteristic time-scales 70, fi. The 
higher cumulants of the distribution, discussed in the next subsection, oscillate and decay 
more rapidly. This is responsible for the inspirals in phase space of the evolution of Gaussians 
plotted by Unruh and Zurek [")], which are calculated in the next subsection and plotted in 
Fig. 6. Of course this is all assuming a non-vanishing ohmic term in the spectral function. 
If the spectral function is purely non-ohmic {i.e., with 70 = 0) then this factor will describe 
the initial state oscillating with a renormalized frequency, but which never decays away. 

The birth factor describes the complicated birth of and settlement into a state of thermal 
equilibrium. This factor is always Gaussian with a covariance matrix given by Eq. (IV.20), 




(IV.21) 
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which involves a convolution of the diffusion matrix with propagators that reflect the natural 
oscillatory decay of the system. This covariance matrix vanishes at the initial time and tends 
at late times to the equilibrium covariance matrix (IV. 21), with the diffusion coefficients 
being their asymptotic constant values given by Eqs. (II.60)-(II.61). This covariance matrix 
is positive definite for all reasonable diffusion constants. Moreover, the anomalous diffusion 
coefficient actually acts as an "anti-diffusion" term that makes axx (and the uncertainty in 
x) free of the log(A/r2^.) divergence, as it will be discussed in the next subsection. 

B. Analysis of the Solutions 

1. Trajectories of the Cumulants 

As we have already mentioned, the Fourier transfom of the reduced Wigner function 
corresponds to its characteristic function, from which the correlation functions for the phase- 
space variables can be easily derived using Eq. (IV.6). The general expressions for the 
cumulants can be obtained straightforwardly from the logarithm of the reduced Wigner 
function in Fourier space as follows: 

oo , n 

E W n^^" = ^^smt, 1^), (IV.22) 

n=l ■ 1=1 

where /c*' denotes the components of the vector k and we used the Einstein summation 
convention for pairs of repeated indices {i.e., it is implicitly understood that a sum Yl'ii=i 
should be preformed over each pair of repeated indices ii). k*^"-* is the n*^ cumulant and 
acts as a tensor of order n contracted with n copies of k. Using the result for Writ, k) from 
Eq. (IV. 19) we have 

OO 71 OO 71 

E iAUit)U^k- = Y: ^-K.JO) (e-*^k)^' - ikV,(t)k, (IV.23) 

n=l ■ 1=1 n=l ' 1=1 

where '^j1^\„(0) are the cumulants associated with the initial distribution. Eq. (IV.23) 
implies 

n . . 

^..^Jt) = «^£.,„(0) [e-^^^Y" + 5^,ar{t). (IV.24) 

We can see that the only cumulant with a non-vanishing asymptotic value, which is a 
consequence of the thermal fluctuations, is the covariance matrix (with n = 2). The closely 
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related second momenta of the distribution are given by 

(qq-)(t) = e-*^(qq-)oe-*« + ^^(t), (IV.25) 

where (■ ■ ■ )o denotes expectation value with respect to the reduced Wigner function at the 
initial time. All other cumulants experience oscillatory decay with time scales of 77,70, nfl 
with n being the order of the cumulant. In particular, the expectation value 

(q)(t) = e-*«"(q)o, (IV.26) 

follows a trajectory plotted in Fig. 6, where one can see that the trajectory of the expectation 
values (x), (p) for any initial distribution inspiral into the origin. 



<P> 

M Qr 




<x> 



FIG. 6: The trajectory of the expectation values (x), {p). 



2. Late Time Uncertainty Function 

As we have seen above, any specific features of the initial distribution decay away (as- 
suming a non-vanishing ohmic term in the spectral function) and at late times the state 
tends generically to a Gaussian with a covariance matrix given by Eq. (IV.21). There- 
fore, from Eq. (IV.25) it follows that at late times (Ax)^ = (cr^)^.^. and (Ap)^ = (cr^)pp. 
From Eq. (IV.21) and using Eqs. (II.32)-(II.33) we can express the position and momentum 
uncertainties at late times as 

(Ax)^ = ^ JfIi, (IV.27) 
(Ap)2 = -70MFI3. (IV.28) 

TT 
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The product of the two uncertainties 

{Ax)\Apy = (^^^ FI1FI3, (IV.29) 

can be expanded for high temperatures as 

{AxnApY =(J^J + i...)T+{... )T°. (IV.30) 

Inspecting the terms in powers of T immediately reveals the high-temperature result of 
classical statistical mechanics for the case of an ohmic spectrum. Integer subohmic terms 
would strictly decrease this amount. One can also see that at weak coupling the uncertainty 
function agrees with the weak coupling approximation for moderate values of the cut-off 
scale, as shown in Fig. 7. 




0.20.40.60.8 1 1.21.4 Qr 

FIG. 7: Late time AxAp for — high temperature, classical statistical mechanics, ■ • ■ weak coupling 
approximation ^ coth - HPZ at A = lO^Jl^, and • • • HPZ at A = 10^ fir- 

Had one naively tried to have finite diffusion coefficients in the limit A ^ 00 subtracting 
by hand the log{A/Qr) term, one would find a violation of the Heisenberg uncertainty 
principle at low temperature and strong coupling (see Fig.8), which renders the theory 
unphysical. Of course this does not happen with the unsubtracted theory, as seen in Fig. 9. 
It is thus clear that the logarithmic dependence on the ultraviolet cut-off that appears in 
the diffusion coefficients is a physically important parameter and not something that can be 
subtracted away. 

From Eqs. (II.40)-(II.41) [or alternatively from Eq. (11.34)] one can see that mentioned 
earlier, only the momentum uncertainty contains a logarithmically divergent cut-off depen- 
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Subtracted Theory 




FIG. 8: Late time AxAp for the subtracted theory. 
A=10-^Qr Unsubtracted Theory 



AxAp 




FIG. 9: Late time AxAp for the unsubtracted theory. 

dence. In contrast, the position uncertainty is much smaller and finite in the limit A — >■ oo 
(this had already been noticed for Gaussian wave-packets in Ref. [j]). In fact, if the Brow- 
nian particle is coupled to a reservoir at arbitrarily low temperature and arbitrarily strong 
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coupling, the uncertainty in position will be arbitrarily small: 

lim limAx = 0, (IV.31) 

70^00 T-»0 

lim lim An = 00. (IV.32) 

Hence, strong coupling to a low temperature reservoir tends to localize the Wigner distribu- 
tion in position. From this result, one would qualitatively expect strong coupling to suppress 
quantum tunneling while high temperature would aid both tunneling and classical escape 
(thermal activation). 



3. Linear Entropy 

In this subsection we investigate the linear entropy [1!)], which can be easily obtained 
from the Wigner distribution as follows: 

5l = 1 - Tr(p2) = l-2n J £qW^{q). (IV.33) 

In Fourier space it becomes 

^L = l-^y" rf2k|W,(k)|2, (IV.34) 
and using the result in Eq. (IV. 19) we finally get 



Sl = 1-^ j d^kWr (0, e-'^k 



At the initial time the linear entropy is that of the initial state, and at late times it tends 
to Sl = 1- (l/2)(detcr5?)-i/2. 

Alternatively, one can express the linear entropy in terms of an integral of the Fourier- 
transformed reduced Wigner function at the initial time by introducing the change of vari- 
ables ko = e '^k. Eq. (IV.35) can then be written as 

Te+*«Crr(t)e+'H^ko 



Sl = l~^J rf'kodet(e+*"")|W.(0,ko)|V'^° 

yl j rfX|W.(0,ko)pAr(o,e+*«o-T(t)e+*^;ko), (IV.36) 



^ 2det^/'[o-r(t)] 



^ In this limit one needs to use the expressions for the overdamping regime, as explained in footnote 3. 
Moreover, one needs to keep A > 70, which means that A also tends to infinity as 7 00. Alter- 
natively, one should use Eqs. (II.40)-(II.41) plus the terms involving negative powers of A which were 
neglected there, so that the restriction A > 70 does not apply. In both cases one gets the result quoted in 
Eqs. (IV.31)-(IV.32). 
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where A^(/x, cr; ko) is a normalized Gaussian distribution for the variable ko with mean fx 
and covariance cr. For small times this integral is similar to that for the initial state, whereas 
for long times the normalized Gaussian distribution becomes increasingly close to a delta 
function. 

For a Gaussian initial state VVr(0,ko) = exp[— ko^cr5-(t)ko — zko(q)o] the integral in 
Eq. (IV. 35) can be explicitly computed: 

2n J 

= I ^ (ly 37) 

2det^/'[e-*Hcroe-*«^ + crr(t)]' ' ^ 

For these Gaussian states, reasonable linear entropy is synonymous with reasonable un- 
certainty functions {i.e., the linear entropy will be positive if and only if the Heisenberg 
uncertainty principle is satisfied). We have already found that the late time uncertainty is 
well behaved. The uncertainty at the initial and intermediate times should not violate the 
Heisenberg uncertainty principle either. As a particular example, let us consider an initial 
state that corresponds to the ground state associated with the bare frequency. At the initial 
time the determinant is completely specified by the determinant of ctq, which equals 1/4 as 
a consequence of the cancellation between the large (Ap)^ factor (of the order of Mfibare 
with i^bare ~ V^) and the small (Ax)^ factor (of order 1/Mf2bare)- However, after a very 
short time, when crxit) starts acquiring non- vanishing values, it will dominate the total 
contribution to (Aa;)^ and the linear entropy will increase suddenly (this kind of behavior 
for the entropy was found in Ref. [20]). At later times it will relax to its thermal equilibrium 
values. 



C. Solutions of the General Master Equation 

In situations more general than the Laurent series spectrum and where the classical equa- 
tion of motion involves fractional calculus {i.e., it becomes an integro-differential equation 
rather than an ODE), the master equation can be of a slightly more general nature, with 
additional time dependence in the frequency and the dissipation coefficients: 



d_ 

Ft 



iy.= (-^^ + M(],^(t)x| + 2r(t)|p + v^,D(t)v,)iy.. (iv.38) 
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The method for finding its solutions will be almost the same, except that solving for the 
characteristic curves of the Fourier-transformed phase-space variables is now less straight- 
forward due to the time dependence of H, and the procedure followed to solve Eq. (IV. 12) 
needs to be generahzed. The system of ODEs that corresponds to 

- = k-H(s), (IV.39) 

can be decoupled through a process of differentiation and substitution that leads to the 
equations for a pair of independent parametric oscillators where solving one effectively solves 
the other up to an integral. The second order ODE satisfied by kp is 

= kp- 2Tkp + (nl - 2f ) kp . (IV.40) 

After solving it, kr^ can be simply obtained as 

k^ = kl\M [ ds' n%s')kp{s'). (IV.41) 
Jo 

Alternatively, one can proceed in the reverse order as follows: 

= k,-(^2T + S j k + , (IV.42) 

Solving for kp first appears to be simpler since everything is then linear in the coefficients. It 
is convenient to factor out the exponential growth, which reduces the ODE to its undamped 
form: 

kp = e^o^s'W^^^ (IV.44) 
= jp+(nl-T^-t)jp. (IV.45) 



If one is able to solve this differential equation, which for simple enough functions could be 
performed with Floquet analysis and variation of parameters, then the characteristic curves 
can be expressed in terms of the following matrix equation: 

k = $(s)ko, (IV.46) 

where is the so-called transition matrix, which has many of the properties of an expo- 
nential and contains the parametric oscillatory behavior of the characteristic curves. If r(t) 
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and flr(t) tend to some asymptotic values at large times, then the behavior of tends 
to that of a damped harmonic oscillator. 

Once the transition matrix is available one can apply exactly the same approach as above 
in order to solve the master equation. One starts by writing the master equation as 

^>V,(t(s), *(s)ko) = -k-S$"(s)D(s)$(s)ko WMs), $(s)ko). (IV.47) 

Next, one solves this linear ODE in s and, reexpressing the result in terms t and k, one can 
finally write the solution of the master equation as 

Wr (t, k) = Wr (0, $"'(t)k) e-ik^'^TWk^ ^jy_4g) 

with 

(Trit) = 2 (^^ [$-^(t)f $^(s)D(s)$(s)$-^(t)rfs^ . (IV.49) 

For ohmic-like coefficients, the qualitative behavior of the solution (IV. 48) is very sim- 
ilar to that of the ohmic solution. The moments of the initial distribution will experience 
oscillatory decay and the solution will tend to a thermal equilibrium state at late times. In 
purely supraohmic regimes where F vanishes, the initial state would remain with renormal- 
ized frequency. In purely subohmic regimes where F grows quickly, the initial state would 
vanish quickly. 



V. INFLUENCE OF A CLASSICAL FORCE 



In this section we consider the case of a classical force F{t) acting on the quantum 
oscillator. This is done by introducing a time-dependent potential —F{t)x: 

L^ = ^M {x^ -n^x"^) +F{t)x. (V.l) 



A. The Master Equation Coefficients 



To derive the form of the master equation we follow the method of Calzetta, Roura, and 
Verdaguer [18] as it is of a very general nature and can be adapted quite straightforwardly. 
In the presence of an external force the system action takes the following form in terms of 
the "center of mass" and "relative" coordinates, X = [x + x')/2 and A = x — x': 



S{x) - S{x') 



M 



Mds 



X{s)K{s) - n''X{s)A{s) + F(s)A(s) 
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(V.2) 



The reduced Wigner function can still be expressed as an average over a stochastic process 
and a distribution for the initial conditions as follows: 

Wr{x,p,t) = /(5(X(t)-x)(5(MX(t)-p),) \ , (V.3) 
[ V^...e-'^^-~'<, (V.4) 



where we used the notation ■ = in the last equation, and X(t) is now a solution to the 
Langevin equation with the external force: 

M + 27o| + fi^) X{t) = F{t) + m- (V.5) 

Differentiating with respect to time reveals the usual dissipation, renormalized harmonic 
potential and the classical force potential —F{t)x: 

^W. = (v^^Uci-Fit)^^Wr-^(^{mHXit)-x)5[MXit)-p))^^ 

After a functional integration by parts, the last term gives rise to the diffusion terms, with 
the following coefficients: 



1 /•* 
^ Jo 



6X{t) 



(V.7) 



Eq. (V.5) can be solved in the same way as in the case without external force to obtain 
X{t). It is simply given by a homogeneous solution that contains all the information on 
the initial conditions, plus a convolution of the external force with the retarded propagator 
associated with the homogeneous part of the equation, plus a convolution of the stochastic 
source with the same retarded propagator. By the linearity of the functional derivative, 
only depends on the retarded propagator, which is the same as in the case of no external 
force. Therefore, the master equation has exactly the same form, except for the addition of 
the classical force in the potential derivative term of the Poisson bracket: 

d 



(V.9) 



dt 

B. Solutions of the Master Equation 

Fourier transforming the phase-space variables, the master equation becomes 

+ k^HVk^ Wr = - (k^Dk + iF{t)kp) Wr ■ (V.IO) 
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This equation can also be solved via the method of characteristic curves, following the same 
method of the previous section: 

^W,(t(s),k(s)) = -{mk + tF{s)kp)Wr{t{s)Ms)), (V.ll) 
^W.(3,e^^ko) = -(k-oe^«D(.)e^^ko + ^F(.)ky^ko)w.(.,e^^ko),(V.12) 

with the solution 

Wr[s,e^^ko)=WAOM)e-^'''<''^"''''''^'^''^^^^^ (V.13) 
which can be finally written as 

Wr (t, k) = Wr (O, e-*^k) e-i'^^'^^W'^"<<'^>^W)^'^ , (V.14) 

where 

(k)^(t) = /" rfs F(s)e(^-*)%. (V.15) 
Jo 

One can see that just as all temperature dependence only appears in the second cumulant, 
or covariance, the external force only affects the first cumulant, or mean. The change in the 
evolution of the mean, given by Eq. (V.15), simply corresponds to adding the response to 
the external force. As with an ordinary driven oscillator, the response is a convolution of 
the retarded propagator, which exhibits the natural oscillatory decay of the system, and the 
driving force. 

VI. SUMMARY OF RESULTS 

Quantum Brownian motion of an oscillator coupled to a thermal reservoir of quantum 
oscillators has been the canonical model for studying the environmental effects on a quantum 
system, even of macroscopc scale, such as quantum dissipation, diffusion, decoherence and 
entanglement. It also provides important information on quantum measurement, such as 
noise, fluctuations, correlations, uncertainty relation and standard quantum limit in meso- 
scopic systems. Many experiments have been carried out for testing these processes. Fifteen 
years ago an exact master equation [12] for the reduced density matrix of the system were 
derived for a general environment of arbitrary spectral density and temperature. There are 
claims that exact solutions have been found [JJ]. In this paper we report on solutions to 
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this equation for a fairly general set of physical conditions and analyze its salient features. 
We expect these solutions to be useful in realistic settings for the analysis of many problems 
which can be described by this model. 

The first question we addressed concerns the initial time divergences. These divergences 
result from the sudden coupling of an initially uncorrelated system and environment. This 
situation is unphysical and can be resolved by switching on the system-environment interac- 
tion smoothly over a time-scale longer than the inverse of the characteristic UV cut-off scale 
of the environment. Having switched on the coupling, the initial frequency of the oscillator, 
as it quickly evolves from the bare to the renormalized frequency, is the only remaining 
initial singular behavior. But this merely implies a transition from the bare initial state to 
a "renormalized" version of that state, which can be taken as the effective initial state for 
the subsequent regular evolution. 

We then solved the master equation coefficients for a general class of bilinear system- 
environment couplings, which includes any combination of ohmic and integer supraohmic 
and subohmic spectral densities with extreme cut-offs. Although it is by no means a com- 
plete class of spectral densities, it is a privileged class in that the classical trajectories are 
determined by ordinary differential equations and not integro-differential or fractional dif- 
ferential equations. In that sense, their dynamics is more classical. 

For these system-environment couplings, we have solved the master equation coefficients 
for all temperature ranges, beyond weak coupling, and at all times after the short initial time- 
scale during which the coupling is switched on. This has enabled us to clarify the validity 
and shortcomings of previously obtained approximations and provide more general results. 
Perhaps the most useful are our simple expressions for the master equation coefficients at 
late time [Eq. (II.60)-(II.61)]. They clearly reveal the existence of logarithmic divergences 
in the limit of infinite UV cut-off not only in the anomalous diffusion coefficient, but also 
in the normal one. These divergences cannot be consistently subtracted and condition the 
possibility of obtaining meaningful results upon the existence of a physically well-motivated 
cut-off. 

For the system-environment couplings under discussion, which all have local dissipation, 
we have also solved the master equation in its entirety. All initial states evolve into a given 
Gaussian state corresponding to thermal equilibrium. We have obtained the covariance 
matrix for this state [Eq. (IV. 21)], and the corresponding uncertainty function (Sec. IV B 2). 
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Interestingly, arbitrarily strong coupling to a zero temperature reservoir will tend to localize 
the position with arbitrary precision. 

For more general systems with non-local dissipation, we have drastically reduced the 
task of solving the master equation to that of solving a one-dimensional classical parametric 
oscillator problem. The master equation solutions are parametrically similar to the case of 
local dissipation. The same kind of terms arise, including a thermal covariance, but it is not 
guaranteed that the system will relax to a thermal state. That of course depends upon the 
behavior of the classical trajectories. 

Finally, we extended the model of the quantum oscillator linearly coupled to a thermal 
reservoir of oscillators by including a classical driving force. This modifies the dynamics by 
driving the mean position and momentum around just as with a classical driven system. In 
this model we found that the force has no effect upon the width of the wave-packet or any 
cumulant other than the mean. These results may be useful for the study of low-temperature 
measurements of forced oscillators, which are relevant for experiments with nanomechanical 
resonators [23, 24]. They also play a crucial role in future schemes for the detection of 
gravitational waves with high-intensity laser interferometers, where the radiation pressure 
effects on the cavity mirrors are important [25, 26]. 
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APPENDIX 



1. Harmonic Number 

The Harmonic Number H(n) is a function similar to a logarithm. Its analytical continu- 
ation to the complex plane is similar too. 

H(n) = X^p neZ+ (A.l) 
fc=i 

H(0) = 0, (A.2) 

jE = lim(H(n) -log(n)), (A.3) 

n— >oo 

R{z) = -fE + Hz + l), zeC (A.4) 

m = (A.5) 

i){z) ^ \nz- — - + ■■■ if |arg(2;)| < TT, (A.6) 
where Eq. (A.6) is the asymptotic expansion for \z\ oo. 



2. Exponential Integrals 

The following definition and properties of the exponential integral are used in this paper: 

/•oo -t 

Ei(x) = -P.V. / dt — if a;>0, (A.7) 

J-x t 

roo -t 

Ei(z) = / dt — if |arg(;z)| < vr, (A.8) 

J z t 

Ei(z) ~ + l + (A.9) 



z \ z z'^ 



where P.V. denotes the Cauchy principal value, corresponding to the singularity at t = of 
the integrand in Eq. (A.7), and Eq. (A.9) is an asymptotic expansion for \z\ oo. 
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